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Abstract 

We report a detailed study of the electronic and structural properties of 
the 39K superconductor MgB2 and of several related systems of the same 
family, namely Mgo.5Alo.5B2, BeB2, CaSi2 and CaBeSi. Our calculations, 
which include zone-center phonon frequencies and transport properties, are 
performed within the local density approximation to the density functional 
theory, using the full-potential linearized augmented plane wave (FLAPW) 
and the norm-conserving pseudopotential methods. Our results indicate es- 
sentially three-dimensional properties for these compounds; however, strongly 
two-dimensional cj-bonding bands contribute significantly at the Fermi level. 
Similarities and differences between MgB2 and BeB2 (whose superconduct- 
ing properties have not been yet investigated) are analyzed in detail. Our 
calculations for Mgo.5Alo.5B2 show that metal substitution cannot be fully 
described in a rigid band model. CaSi2 is studied as a function of pressure, 
and Be substitution in the Si planes leads to a stable compound similar in 
many aspects to diborides. 

PACS numbers: 74.25.jb, 74.70.-b, 74.10.+V, 71.20.-b 
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I. INTRODUCTION 



The recent discovery of superconductivity at ~ 39K in the simple intermetalhc com- 
pound MgB2 is particularly interesting for many reasons. In the first instance, this critical 
temperature is by far the highest if we exclude oxides and Ceo based materials. Furthermore, 
B isotope effect suggests that MgB2 is a BCS phonon-mediated superconductor, with Tc 
above the commonly accepted limits for phonon-assisted superconductivity. Another reason 
of interest is that the AlB2-type structure of MgB2 (which can be viewed as an intercalated 
graphite structure with full occupation of interstitial sites centered in hexagonal prisms 
made of B atoms) is shared by a large class of compounds (more than one hundred), where 
the B site can be occupied by Si or Ge, but also by Ni, Ga, Cu, Ag, Au, Zn and Cd, and 
Mg can be substituted by divalent or trivalent sp metals, transition elements or rare earths. 
Such a variety is therefore particularly suggestive of a systematic study of diborides and 
structurarly related compounds. To confirm this interest, a recent study 0] showed that 
superconductivity is supressed by addition of more than about 30% of Al atoms. Pressure 
effects also lead to a depression of superconductivity. 

Early studies of the electronic structure of diborides by the OPW and tight-binding 
methods have shown similarity with graphite bands. A more systematic study of AlB2-type 
compounds by the full-potential LAPW method showed in simple metal diborides the 
presence of interlayer states, similar to those previously found in graphite [|lO[|lll, in addition 
to graphite-like a and vr bands. Transition metal diborides, on the other hand, show a 
more complex band structure due to the metal d states. A very recent electronic structure 
calculation [|12| on MgB2 shows that metallicity is due to B states. Phonon frequencies 
and electron-phonon couplings (at the F point) lead these authors to suggest phonon-based 
superconductivity. Further, new calculations |jl3| on related compounds (LiB2 and graphitic- 
like B) addressed the peculiar contributions to conductivity coming from the interplay of 
2D and 3D Fermi surface features. An alternative mechanism has been proposed by Hirsch 



\TM, based on dressed hole carriers in the B planes, in the presence of almost filled bands. 
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Recently, superconductivity at Tc ~ 14K was observed in CaSi2 [|T5|. In this compound, 



high pressure has been effectively used to tune both structural and physical properties 
T^,[T6|. In particular, Bordet and co-workers [|T^ identified a new structural phase transition 



for CaSi2 from trigonal to hexagonal AlB2-type structure, with superconductivity at Tc = 
14 K, one of the highest transition temperatures known for a silicon based compound. The 
occurrence of this structural phase transition indicates that Si shifts from an sp^ to an sp"^- 
like hybridization, and experiment points to a clear enhancement of Tc with high pressure 
Tsj . Many electronic structure studies have been carried out for the different phases [[17| -[T9 



of disilicides. CaSi2 in the AIB2 structure was studied by Kusakabe et al. |2^, but a final 
answer on the stability of this phase under pressure is still lacking. 

As previously pointed out, the AIB2 family of compounds provides a wide playground 
for the search of new superconductors, and we intend to provide here a contribution, based 
on a detailed study of electronic and structural properties together with phonon spectra, 
for several diborides and for CaSi2, in the trigonal and hexagonal structures. The questions 
that we want to address are the following: (i) how do the electronic properties of MX2 
depend on the atomic species M; namely, can rigid band schemes provide the correct trend, 
or is the metal effect mostly mediated by lattice parameter changes (chemical pressure), (ii) 



Ref. shows the presence of both B vr Fermi surfaces (FS) giving rise to 3-dimensional 
states, and of small 2-dimensional B a-bonding pockets. It is of course interesting to know 
which of the FS is mostly involved in superconductivity. The trends observed in the band 
structures, together with experimental studies, will contribute to elucidate this point. In 
the same way, changes in the Fermi level density of states with chemical composition will 
possibly correlate with Tc values, (iii) Pure CaSi2 shows a trigonal instability, joined by a 
buckling of Si2 planes, which is apparently ||T^ strongly reduced under high pressure, or by 
Be substituting some of the Si sites. The presence, or the proximity, of structural instabilities 
(investigated through phonon spectra), will hopefully provide hints to further experimental 
search for new materials. 

In this paper, we present calculations for MgB2, Mgo.5Alo.5B2, BeB2, CaSi2 and 



CaBexSi2-x {x = based on the local density approximation (LDA) to the density func- 
tional theory. We use both the Full Linearized Augmented Plane Wave method, and the 
norm-conserving pseudopotential method. In Sect. |l|we give some computational detail; in 
Sect. |T| and ^ we discuss the electronic properties for all the compounds considered; we 
finally draw our conclusion in Sect. |V]. 



II. COMPUTATIONAL AND STRUCTURAL DETAILS 

Our calculations were performed using the "all-electron" full potential linearized aug- 
mented plane waves (FLAPW) ||2l| method, in the local density approximation (LDA) to 
the density functional theory. In the interstitial region we used plane waves with wave vector 
up to Kmax = 3.9 a.u. and 4.1 a.u. for disihcides and diborides, respectively (we used the 
larger K^ax = 4.4 a.u. for BeB2 because of the smaller spheres imposed by the smaller 
lattice parameters). Inside the muffin tin spheres, we used an angular momentum expansion 
up to Imax = 6 for the potential and charge density and Imax = 8 for the wave functions. 
The Brillouin zone sampling was performed using the special fc-points technique according 
to the Monkhorst-Pack scheme and also the linear tetrahedron method with up to 120 



A;-points in the irreducible Brillouin zone. We used muffin tin radii Rmt = 2.2, 1.8, 1.59 
and 2.1 a.u., for Ca, Si, B, and the remaining elements respectively (we used Rbc = 2 a.u. 
and Rb = 1-5 a.u. in BeB2). 

The crystallographic structure of AIB2 (C32) is hexagonal, space group PQmmm with 
graphitic B planes, and 12-fold B-coordinated Al atoms sitting at the center of B hexagonal 
prisms. The unit cell contains one formula unit, with Al atom at the origin of the coordi- 
nates and two B atoms in the positions (di = (a/2, a/2-\/3, c/2), d2 = — di). The trigonal 
distortion in CaSi2 corresponds to a buckling of Si planes leading to a P3ml space group 
symmetry), with one internal structural parameter z (the internal coordinate of B atoms). 

To calculate plasma frequencies, Fermi velocities. Hall coefficients, and the densities of 
states (DOS) near Ep, we have performed a spline fit of the ab-initio energy bands. 
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The fitted bands have then been computed over 128^ points in the reciprocal space unit cell. 
The linear tetrahedron method has been used to compute DOS and the other quantities. 

Phonon frequencies were computed using the ABINIT code by means of Linear 
Response Theory ||25|J2^ (LRT) approach in the framework of pseudopotentials plane wave 
method |^^. Troullier and Martins soft norm-conserving pseudopotentials were used for 



all the elements considered, Mg and Ca p semicore states were not included in the valence, 
their contribution to the exchange and correlation potential was accounted by the non- 



linear core corrections [^]. The local density approximation to the exchange-correlation 
functional is used. Plane wave cutoffs of 44 Ry and 30 Ry were found sufficient to converge 
structural and electronic properties of the diborides and silicides respectively. Brillouin zone 



integration was performed using a 20 x 20 x 20 Monkhorst-Pack p2[ mesh for the structural 



optimization and a 12 x 12 x 12 mesh was used in the LRT selfconsistent calculation of 
phonon frequencies. To improve fc-point sampling convergence a gaussian smearing with an 
electronic temperature of 0.005 Ha was used. The structure used in the LRT calculations 
was obtained optimizing the lattice structure up to a residual stress of 1 MPa. 

III. DIBORIDES 
A. MgB2 and BeB2. 

We shall start our discussion of electronic properties of diborides by presenting in Fig. |l| 
the energy bands of the superconducting compound MgB2, which will be the reference for 
all our further studies (the M — T — K lines are in the basal plane, while the L — A — H lines 
are the corresponding ones on the top plane at kz = tt/c). We used the experimental lattice 
parameters a = 3.083 A, c = 3.52 A. Our results, which are in excellent agreement with those 
of Kortus et. al. [jl2[, show strong similarities with the band structures of most simple-metal 



diborides 0, computed by the FLAPW method. They show strongly bonded sp"^ hybrids 
laying in the horizontal hexagonal planes, forming the three lowest (cr-bonding) bands, as 
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the main structure in the valence region. The corresponding antibonding combinations are 
located around 6 eV above the Fermi level {Ep). 

The bonding along the vertical direction is provided by the tt bands (at ~ -3 eV and ~ 
2 at r), forming the double bell-shaped structures. Similarly to AlB2-type diborides, and 
unlike graphite, we notice a large dispersion of the vr bands along the direction {T — A). 
This can be understood on a tight-binding approach: the phases of the Bloch functions 
lead the B pz orbitals on adjacent layers to be antibonding at F, and bonding at A. The 
downwards dispersion of the vr-bonding band (from —3 to —7 eV along T — A) is joined 
by an opposite upwards dispersion of an empty band, (at ^ 2 eV at the F point). The 
analysis of the wavefunctions for this empty band reveals that this band has an interstitial 
character, and is very similar to those found in graphite in diborides 0], in CaGa2 



0, but also in CaSi2. A general description of the bonding can be given by the partial 
density of states (PDOS), shown in Fig. The B PDOS, in particular, shows bonding 
and antibonding structures for the s and p states. The Pz states give rise to a rather wide 
structure responsible for less than half of the DOS at Ep, which therefore results having 
predominantly B p — a character. 

In the search of similar materials having the desired superconducting properties, BeB2 
is the first natural candidate, as the band filling is expected to be pretty similar, and lighter 
Be atoms may help providing larger phonon frequencies while keeping similar electronic 
properties. The first experimental report on BeB2 gives the average lattice constants. 



These experimental parameters have been used for all the previous calculations of the BeB2 
0J^. We have optimized the lattice constants, imposing the AIB2 structure, obtaining a = 
2.87A, c = 2.85A, smaller than in MgB2 and similar to the values estimated by the average 
experimental values. In particular, we notice that while a reduces only by about 7.4%, c 
reduces by about 24% when Mg is substituted by Be, a result which can be understood 
pretty easily. On one hand, in fact, the optimization of cr-bonds prevents a too drastic 
variation of the a parameter; on the other hand, c can change more easily as the vertical 
bonding, provided by the vr bands, has a large contribution from the metal orbitals. The 



band structure of BeB2, shown in Fig 0, are pretty similar to those computed by the FLAPW 
method |Q| at the experimentally estimated lattice parameters. We can remark the strong 
similarities with the MgB2 compound, which includes the presence of a-bonding hole pockets 
along r — A, but also a few relevant differences. First of all, the shorter lattice parameters 
lead to wider valence bands, and also to more dispersed a bands, in particular the cylindrical 
hole pocket along T — A, which might be relevant for superconductivity. Also, the different 
energy location of the metal s-free electron band (lower in MgB2), and the different shape 
of the vr-bonding band (related by the different c values) causes a different occupation of 
this band, especially along F — M. 

To further investigate these two materials, and their relation to isostructural diborides, 
we studied the zone center phonon frequencies of MgB2, BeB2, and AIB2, using the linear 
response approach, with the pseudopotential method. The calculations have been carried 
out at the theoretical lattice constants (also listed in Table D, as obtained using the same 
method. The two lowest independent frequencies correspond to a motion of the intercalated 
metal atoms relative to the rigid B networks: one vertical mode {A2u), and one doubly 
degenerate in-plane mode, (Eiu). The two higher frequencies represent the internal motion 
of the B network itself: again one vertical mode (-Big), now representing a trigonal distortion 
producing a buckling of the planes similar to that found in CaSi2, and one degenerate {E2g) 
mode giving an in-plane stretching of B-B bonds. Our calculations for MgB2 are in good 
with other results |ll2| , p0| . The comparison between the MgB2 and BeB2 shows a much 
higher frequency in BeB2. This result cannot be explained through the different atomic 
masses, as it involves only B motion. It is therefore related to the different coupling of the 
electrons with the nuclear system. Surprisingly, two of the remaining frequencies are lower 
in BeB2, namely the A2U mode which involves the motion of lighter Be atoms against the 
B network and the Big mode. This indicates a strong dependence of the phonon spectra 
of these materials on the electronic structure details. The lowest Eiu mode has comparable 
frequencies. If we now go to AIB2, we see again that the E2g mode is much higher than in 
MgB2. Therefore, this mode has his minimum value in the superconducting material. The 
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modes involving the relative motion of the B2 planes as a whole have frequencies comparable 
to the other materials; the Big, on the other hand, is smaller than in MgB2 and BeB2, which 
is probably a consequence of the larger band filling (see the discussion on CaSi2 below). 

We have computed, within a rigid band scheme and using the scheme described, e.g., in 
Ref. the Fermi velocities, plasma frequencies (along the principal axes of the crystal). 



and the independent components of the Hall tensor for BeB2 and MgB2. Because of hexag- 
onal symmetry, the non-zero components will be Rxyz, corresponding to the magnetic field 
along the z axis and transport in plane, and Rzxy = Ryzx corresponding to in-plane mag- 
netic field. The results, plotted in Fig. ^ as a function of the hole or electron-type doping, 
show relatively smooth variations over a wide range of doping, and are nearly coinciding, at 
zero doping, with the corresponding calculations by Kortus et al. |12[. After comparing the 
two materials, the following remarks can be done: (i) the density of states at Ep is smaller 
in BeB2. (ii) MgB2 is nearly isotropic in terms of plasma frequencies and Fermi velocities, 
while BeB2 is expected to show important anisotropy in resistivity, and, in particular, the 
in-plane conductivity is expected to be smaller. A band-by-band decomposition explains 
these results in term of very similar (and very anisotropic) contributions in the two com- 
pounds from (J states, and quite different contributions from the vr bands, nearly isotropic 
in MgB2 and favoring conductivity along the c-axis in BeB2. These results point out that 
despite structural similarity, the electronic properties of AlB2-type structures differ substan- 
tially from graphite and intercalation graphite compounds, (iii) The Hall coefficients are 
positive (hole-like) when the magnetic field is along the z axis (transport in-plane), while 
is negative for MgB2 and almost zero but negative for BeB2, when the magnetic field is 
in-plane. A detailed analysis shows again comparable contributions in the two compounds 
from a bands, while the sign of Rzxy and Ryzx comes from a balance of the vr bands. The 



Hall coefficient of MgB2 has been recently measured by Kang et al. |31[]. Rh, which will 
correspond in this case to an average of the tensor components, turns out to be positive 
and to decrease with temperature. At T = lOOK they give Rh = 4.1 x lO^^^m'^/C, and 
an extrapolation of their results at T — > 0, gives Rh ~ 6.5 x lO^^^m'^/C. If we average, at 
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zero doping, our tensor components, the positive contributions overcome the negative ones, 
resulting in an average value Rh ~ 2 x lO^^^m^/C. This value is thus of the correct order, 
but single crystal measurements are probably necessary to have a significant comparison, 
avoiding an average of quantities with different sign. It is worth noticing that, if the sign of 
carriers has a role on superconductivity [Q, BeB2 also has a positive Hall coefficient. 

B. Al doping in MgB2. 

A recent experimental report by Slusky et al. 0] shows that Al doping destroys bulk 
superconductivity when the Al content x is greater than ^ 0.3. Stimulated by these results, 
we studied the chemical substitution of Mg with Al. The substitution corresponding to 
a 50% concentration of Al (x = 0.5) leads to a reduction of the cell parameters, from 
a = 3.083 A, c = 3.521 A Q for MgBs, to a = 3.047 A, c = 3.366 A ^ for Mgo.5Alo.5B2. 
These differences again indicate hardly compressible B-B a bonds, and interlayer distances 
quite sensitive to the intercalated cation. The study of Slusky et al. shows consistent 
results, and indicates, furthermore, that after a two-phases region for Al substitutions from 
10 to 25%, c collapses and bulk superconductivity disappears. We first studied this problem 
in a rigid band approach, bringing MgB2 to the same lattice parameters of Alo.3Mgo.7B2 
0]. The only apparent change at that electron concentration is the complete filling of the vr 
bonding band at the M point. The a bonding bands need a larger electron addition to be 
completely filled in a rigid band model. 

To investigate further this problem, we studied the system with 50% Al concentration. 
We used in our calculations the experimental values, and in a first approach simulated 
the disorder (certainly present in the real compound) using an orthorhombic supercell. The 
corresponding unit cell contains two formula units, and is described by the two orthogonal in 
plane basis vectors ai = (a, 0, 0) and a2 = (0, av^, 0). We started by calculating the energy 
bands using the MgB2 lattice constants, and then we used the experimental Mgo.5Alo.5B2 
constants. The results of these calculations are shown in Fig. |^ together with the bands of 
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MgB2, folded into the Brillouin zone (BZ) of the orthorhombic structure. The folding is such 
that the F point of the orthorhombic BZ corresponds to the F and M points of the hexagonal 
BZ, while Z corresponds to both the hexagonal A and L. Therefore, the T — Z line collects 
the T — A and the M — L hexagonal lines, and we notice immediately the three tt bands (two 
bonding and one anti-bonding, coming from folding). The F — X line contains the (folded) 
hexagonal F — M line. If we first look at the frozen structure calculation Fig. ^(c), we see 
quite similar band dispersions with the obvious splittings, related to the symmetry lowering 
induced by the different cations, and a different band filling. We notice, however, a lowering 
of the folded tt antibonding band, related to the different cations: the vr band is in fact 
pushed downwards by the interaction with the free-electron metal-s band, which is at lower 
energy in the compounds with Al (this band is partially filled in pure, non superconducting, 
AIB2 [§]). As a consequence, the extra electrons coming from Al do not fill the a-bonding 
band completely, leaving hole pockets around Z. As we relax the structure, the a bands 
increase their width due to the smaller a value, and the situation near Ep changes only 
because of the occupation of the free-electron band around F. These results suggest that 
although the rigid band scheme can be considered roughly correct, the dependence upon 
doping of the physical properties and of Fermi surface states relevant for superconductivity, 
may be crucially dependent on the actual substitutional atom. 



IV. CaSi2 AND RELATED COMPOUNDS 

When prepared at ambient pressure CaSi2 has a rhombohedral structure (space group 
R3m, a = 3.85 A, c = 30.62 A) and is a semimetal not superconducting down 
to 30 mK 1^^. By annealing under pressure - typical conditions are 8 GPa and 800 K - a 
tetragonal superconducting [a — ThSi2 — type) phase appears with a 3-dimensional network 
P^l of Si, showing superconducivity with = 1.58 K. At pressures between ^ 7 GPa and 
9.5 GPa, rhombohedral structure samples undergo a transition to a trigonal phase (P3ml, 
ax ~ 3.78 A, ct ~ cr/6 ~ 4.59 A) with Si at 2d positions, and z = 0.4, where z is the 
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coordinate of the B atoms in units of c (2; = 0.5 corresponds to the perfect AIB2 structure). 
A second phase transition occurs at 16 GPa. This new phase may be described as [^] an 
AlB2-hke structure with z shghtly smaller than the ideal value 1/2. The AlB2-type structure 



was previously observed in rare-earth silicides ||36|, while in CaSi2 it is apparently stabilized 
by doping the Si planes with Be 

In Fig. P we show the energy bands of CaSi2 in the AIB2 structure (dashed lines), 
superimposed with those of the trigonal structure (full lines), along the main symmetry lines 
of the Brillouin zone. Since this polymorph can be obtained only under pressure, we do not 
minimize the lattice parameters, but rather use the experimental values a = 3.7077 A and 
c = 4.0277 A, corresponding to P=19.3 GPa, and optimize the internal parameter z. The 
bands in the AIB2 structure show strong similarities with those of simple metal diborides, and 
of CaGa2 @|. The major difference is provided by Ca d states, which are mostly found in the 
conduction region, above 4 eV; however, there is a p-d mixing between Si and Ca, especially 
concerning Si vr and Ca d^i orbitals, also mixing with the free-electron states. Unlike in 
MgB2, the Fermi level of CaSi2 cuts only the antibonding part of the vr bands, while the 
bonding part is fully occupied. This of course implies a reduced contribution of vr states 
to the stabilization of an sp"^ environment, consistent with the sp^-like trigonal distortion. 
We will come back to this point later. The energy bands in the trigonal structure show 
that, a part from a small rigid shift, the a bands change very little with respect to the AIB2 
structure. The vr bands, on the other hand, show anticrossings which in the antibonding case 
are located right at the Fermi level. This is most clearly indicated in the density of states 
of CaSi2, reported in Fig. |^. The dip at Ep in the trigonal structure provides a textbook 
explanation for the lattice distortion. 

In order to understand the stability of AlB2-type polymorph, we study the total energy 
variation as a function of z, relative to the ideal z = 0.5 case, for two experimental structures 
p!5| , and report it in Fig. ^(a). There is a good agreement between theory and experiment 
(within 2.5%), and a non-negligible (~ 30 mRy per formula unit) stabilization energy. We 
notice, however, that the stabilization energy decreases with pressure. As we fix a and c to 
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the values corresponding to the larger pressure of 19.3 GPa |T^, we see a larger discrepancy 
relative to the experimental equilibrium value of z [z^xp ~ 0.448). This indicates that some 
different modification might be going on in the experimental samples. In order to discuss 
the link between the stability of CaSi2 in the AIB2 structure and the vr— band filling, without 
modifying the Si2 planes themselves, we vary the atomic number of the cation, which now 
represents a virtual atom ranging from Ca to K (corresponding to Z = 20 + x). Our 
results, reported in Fig H(b), show that while x = —0.3 leaves the total energy curve almost 
unchanged, x = —0.7 or larger give a nearly vanishing trigonal stabilization energy. This 
is consistent with our previous conclusion that a filling of 7r-antibonding bands destabilizes 
the AIB2 structure. 

This last conclusion is supported by the existence, in the hexagonal structure, of 



CaBexSi2-x7 with x = 0.75, at a = 3.94 A and c = 4.38 A [^. In fact. Be doping re- 
moves electrons from the Si2 planes, and brings the formal filling of the above vr-antibonding 
bands to the level found in diborides. To investigate these effects, we studied the (artifi- 
cially ordered, but existing) CaBexSi2-x system for x = 1, not far from the experimental 
stoichiometry. The corresponding bands, shown in Fig. |^, have been calculated at the exper- 
imental lattice parameters ordering Be and Si atoms as, e.g., in hexagonal BN. In this way, 
the Brillouin zone is the same as in the undoped compound, and the interpretation is easier. 
The CaBeSi band structure shown are overall similar to that of the undoped compound. 
The symmetry lowering induced by substitution produces the large splitting of degeneracies 
at the K and H points. The major difference is of course related to the position of the Fermi 
level, which is almost exactly coinciding with that found in the superconducting MgB2. The 
reduced bandwidth for the a bonding bands may be associated with the smaller size of Be 
orbitals. The electron depletion from the antibonding vr-bands removes the trigonal insta- 
bility, as demonstrated by total energy calculations. In fact, while the F- point phonons of 
CaSi2 in the AIB2 structure show the presence of an instable phonon, in agreement with 
FLAPW total energy calculations, in CaBeSi the phonons are stable, as shown in Table |. 
Table H shows that this system has lower phonon frequencies than the diborides, as is normal 
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given the heavier atomic species involved. If we now immagine to vary x in CaBexSi2-x, 
we will therefore have that, at some value of x, the Big mode, corresponding to the trigo- 
nal distortion, will be close to an instability. This situation may thus offer an interesting 
playground in the search for superconducting materials. 

The Fermi level of CaBeSi cuts the a bands, leaving hole pockets quite similar to those of 
MgB2. While it is difficult to say, on the basis of our electronic structure calculations only, 
whether the similarities with MgB2 band structure can lead to superconductivity, it may be 
useful to provide hints on the stability of this compound, which may help experimentalists in 
the search for new compounds. The PDOS in Fig. |10] show a strong hybridization between 
Si and Be states, in the valence region. In other words, substitutions on the Si site can 
be possible without disrupting the sp"^ network if the one-electron energies of the impurity 
atoms are comparable with those of Si, and light atoms should be preferred in order to keep 
high values of phonon frequencies. We may suggest B, which however corresponds to a hole 
doping smaller than Be. 

V. CONCLUSIONS 

We have performed electronic structure calculations for MgB2, BeB2, Mgo.5Alo.5B2, CaSi2 
in the trigonal and in the AlB2-type polymorphs at the high pressure experimental structural 
parameters, and for the CaBexSi2_x system {x = 1), experimentally stabilized for x = 0.75. 
We have calculated band structures, PDOS, transport properties and phonon frequencies, 
using the FLAPW and the norm-conserving ab-initio pseudopotential methods, within the 
local density approximation. The following conclusions can be drawn on the basis of our 
calculations: as compared to MgB2, BeB2 has a similar filling of a bands, but differs in terms 
of TT states near Ep. The DOS at Ep is lower, and some phonon frequencies are substantialy 
higher while others are comparable or lower. It would be interesting to further investigate 
this compound. Al substitution cannot be fully simulated by a rigid band model, as the 
71 antibonding states appear to be quite sensitive to Al substitution, and absorb much of 
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the added electrons. Transport properties show a substantial two-dimensionality of these 
compounds, due to tt states. 

We studied CaSi2 in its trigonal structure, and show that the distortion can be removed 
only by subtracting (in a virtual crystal-like approximation) at least 0.7 electrons. In fact, 
CaBeSi shows no tendency towards trigonal distortion, and electron states near Ep of similar 
bonding nature to those of MgB2. 
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TABLES 

TABLE I. Theoretical structural parameters (in A) and phonon frequencies at F (in cm~^) 
for the AlB2type systems. Level degenaracies are shown in parenthesis. 





MgB2 


BeB2 


A1B2 


CaBeSi 


a 


3.045 


2.903 


2.976 


3.914 


c 


3.480 


2.853 


3.237 


4.488 


(2) E^^ 


349 


327 


306 


189 


(1) A2u 


426 


213 


437 


247 


(2) E^g 


575 


893 


996 


461 


(1) Bu, 


708 


642 


515 


346 
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FIGURES 

FIG. 1. Band structure of MgB2. All the energies are referred to the Fermi level, taken as zero. 



FIG. 2. Total and partial density of states (PDOS) for MgB2. 

FIG. 3. Band structure of BeB2. 

FIG. 4. Transport properties of MgB2 and BeB2 as a function of doping in a rigid band scheme. 
Full and dashed lines correspond to BeB2 and MgB2 respectively. Circles (squares) indicate the 
y (z) principal values of the plasma frequency 17^ (in eV), of the mean squared Fermi velocity vp 
(in 10'' cm/sec) and of Rxyz (in 10^-'^'^ m'^/C, notice that Rzxy = Ryzx)- DOS are in states/eV-cell. 

FIG. 5. Band structure of MgB2 (panel (a)) and Mgo.5Alo,5B2 at the experimental (panel (b)) 
and frozen MgB2 (panel(c)) lattice parameters. 

FIG. 6. Band structure of CaSi2 in the ideal AIB2 (dashed lines) and trigonal (full lines) 
structures . 

FIG. 7. Total and partial density of states (PDOS) for CaSi2 in the AIB2 (dashed lines) and 
in the trigonal structure (full lines). 

FIG. 8. Upper pannel: total energy of CaSi2 as a function of the internal parameter z, for two 
different values of pressure, at the corresponding experimental lattice constants ( a = 3.7554 A, 
c = 4.3951 A for P = 15.0 GPa and a = 3.7668 A, c = 4.4752 A for P = 12.8 GPa, after p|). 
The experimental values of z are indicated by arrows. Lower pannel: same as above, at the larger 
experimental pressure (19.3 GPa), with Ca substituted by a virtual cation having nuclear charge 
Z = 20 + X. The experimental lattice parameters are a = 3.708A and c = 4.028A, z = 0.448). The 
zero of energy corresponds to the AIB2 value z = 0.5. 
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FIG. 9. Band structure of CaBexSi2-x {x = 1) in the hexagonal Brillouin zone, at experimental 
lattice parameters. 

FIG. 10. Total and partial density of states (PDOS) for CaBexSi2-x {x = 1) in the AIB2 

structures (see text). 
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Fig. 1 G. Satta et al. 





Fig. 3 G. Satta et al. 




4 G. Satta et al. 



n (electrons) 



(a) (b) (c) 




z rx sz rx sz rx s 

Fig. 5 G. Satta et al. 
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